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Abstract 

We consider a model of oscillator with nonpolynomial 
interaction admitting exact solutions both for energy 
eigenvalues in terms of zeros of Bessel functions considered 
as * unctions of the continuous index, and for the 
ccrresponding eigenstates in terms of Lommel polynomials 


Let us consider the following Hamiltonian, 


+ . f + ■ + r f ^ / 4 \ 

hi — u Cl ■+■ X 4 CL Tci u. + 1 J + O. + / J Cl > (1) 

Here cl and cl are usual boson annihilation and creation operators, u> 
and X are positive real parameters (the general i zat i on to complex 
coupling constant X does not lead to any new result, since the phase 

L i p 

of X is trivially eliminated by the canonical transformation a -* ae 
preserving the energy spectrum)- If the mean number of quanta is 
close to zero, then (1) turns into the Hamiltonian of usual forced 
oscillator- In the opposite quasi c 1 assi cal regime of large mean 

+ i ..*2 1 ui 

number of excitations N = <a a> » 1 the substitution a ^ AT e 


leads tc the energ y-i ndependent interaction Hamiltonian 


H - Xcos'p , (2) 

W*t 

which is in fact exact, since the expression inside the figure brac- 
kets is nothing but the Sussk i nd-Gl ogower cosine phase operator 111 
which properties were discussed in detail in the known review by 
Car^uthe^s and Nieto "21- 

Expanding the energy eigenstate | E> over the Fock states 


o* f 
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|£> = l< Jn > 


and taking into account the known matrix elements of operators a 
a ore can easily reduce the stationary Schroedinger equation to 
following set of coupled linear algebraic equations, 

Ec = Xc , 

o i 


Ec = u»i + X f c + c 1 , n 2: 1 
n v n-i n+i^ 

It is convenient to introduce dimensionless variables 

£ = \/<jj , /j = £'/o>* (5) 

Then normalized energy /j is determined from the equation $( 2 ,p)= 0 , 
where function f* is the char acter i st i c determinant of system (4): 
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Expanding this determinant over the elements of the first row one can 
easily obtain the following recurrence relation, 

$(£■, fu) = — /j— 1 ) — fj— 2) • (7) 

Introducing new function 

F ( 2 , fj ) = 2 ^5* ( £ , p ) ( 8 ) 

one can rewrite (7) as follows, 

F<*, + F<£, ,u- 2) = -^H-F(£, p-l>. (9) 

But this is the well known relation for Bessel functions C3,43. 
Consequently, the energy levels are determined by zeros of Bessel 
functions in accordance with the equation 

( 2*> = Q. (10) 

For small values of parameter £ the well-known power series expansion 
of the Bessel function leads to the equation 


UL 1 

.ir 

)?i \ T Cm-jjj 


For £ 0 the solutions of this equation with respect to fj are 

determined by the poles of gamma-function- Evidently, they reproduce 
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approximately the equidistant harmonic oscillator spectrum: ~ n 


n = 0,1,2,... Since all poles of gamma-function are simple, with the 

residues (-l) r ’/n!, the correction to the n-th energy 

_ , 2<n+i> 

order of £ : 


level has the 


2 

— -O 


2<ri + l> . 

jj = n — 2 / n 

r* 


1 - 


< 12 ) 


Note that all corrections are negative- 

For large values of the coupling constant we can use the known 

asymptotic formula 


J „ (2.2) * (n2> 


- 1/2 


( 13) 


cos (22 + n^j/2 + rr/4)- 

Then for j/jj « | 2 | the spectrum is equidistant again, but with the 
twice distance between the neighbouring energy levels: 

u 2 s 1/2 + 2 n -4 z/n + Oiz 4 ) - (1^) 

r» 

Here n is an arbitrary integer having the same order of magnitude as 
the 1 arqe parameter 2 - Note that energy values depend on the coupling 
constant in a specific almost periodic manner: 


jj ( 2 ) 


u (2 + n/2 ) - 
r»*M 


( 15) 


Now let us look again at eq- (4)- Comparing it with different 
recurrence relations for special functions given in C*_>,4D, one can 
recocnize that it is nothing but the equation for Lommel ' s polyno- 
mials (which are in fact polynomials with respect to l/2> 


R <T2 1 + R = 

m.i-* r» — i t l- 


k (l-> 


n) 


R 

ri,V 


( 16) 


Consequent 1 y , 


c <Ai> =N(fj,z) R C-2-2.') = 

r, 
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Z C-/J> n " l Cn-lJ>! Ten- 1 -pJ> 
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l \Cn-2L0 \TC l-pJ> 

1=0 


coefficients = c </J> /W(/ j,*> are as follows. 


t o - 1 ’ 


= *i/a, - 1 

A 2 


Taking into account (12) we have, e-g-, for the ground state 


:<o> _ 1 

c. . - 1 9 
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'( o> 


= - 2 , 


'iO> 


= 2 


( 17) 


w^ere ?V(p,2) is the normalizing factor- For example, the first three 


(18) 


(19) 


A '2 

In conclusion let us discuss the correspondence between the 
quantum problem under study and its classical counterpart described 
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: o the energy - phase canonical variables 


E ~ ^C’p 2 + * <P - arccos 2 ] 

with Hamiltonian 


( 20 ) 


H = E + Xc osf . (21) 
Since this Hamiltonian depends linearly on the energy variable E, the 
canonical equations of motion 

3E& t = fiHs&p. dips' dt = -&HS6E (22) 
can be + ound without difficulty for an arbitrary "potential" / <p) : 


pC tJ> = -t * EC' O = E^ + /C"0-> - /C*0- (23) 

However, in the quantum case just the "potential” cos*> seems 
di sti ngui shed- Fo^ example, if one takes instead of (2) the interac- 
tion Ham: 1 tori an 

tm = x [ cQs ^]t (24> 

-v 

then instead of (6> and <7> one gets (fj = /j - 2^) 
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S(jb, = — /j$ (£, ^-1) + £ 2 ( 1 ) 5 r (js , /j- 3) + £ 4 f ( s , jj-4) (26) 

with unknown solution. 

Although the physical meaning of the quantum model with 
Hamiltonian <1) is not clear at the moment (its "nearest neighbour” 
H = E?/2 + Xcos*> describes the Jcsepbson junction), we hope that due 
to its beauty it will find applications in future. 

The authors thank V. P- Karassi ov for valuable discussions- 
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